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Neutral B Meson Mixing and Heavy-Light Decay Constants from 
Quenched Lattice QCD * 

Laurent Lellouch a t a nd C.-J. David Lin b * (UKQCD Collaboration) 
a Theory Division,CERN,CH-1211 Geneva 23, Switzerland 

b Department of Physics & Astronomy, The University of Edinburgh, Edinburgh EH9 3JZ, Scotland 

We present high-statistics results for neutral B-meson mixing and heavy-light-meson leptonic decays in the 
quenched approximation from tadpole-improved clover actions at /3 = 6.0 and j3 — 6.2. We consider quantities 
such as -Bs d(a) , fn d , s s, /s d(s) and the full AB = 2 matrix elements as well as the corresponding S'l7(3)-breaking 
ratios. These quantities are important for determining the CKM matrix element |Vtd|- 



1. INTRODUCTION 



B ( j oscillations allows a clean 



The study of B° d 
extraction of the poorly known CKM matrix ele- 
ment \Vtd\- However, the accuracy of this deter- 
mination is currently limited by the theoretical 
uncertainy in the calculation of the matrix ele- 
ment, 

M bd = (B°\Or= 2 \B° d ) 

= (B° d \h P (l- l5 )dbY(l- l5 )d\B° d ) , 

which is related to the mass difference of the two 
mass eigenstates of the B d — B® system, 
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where Gf is the Fermi constant, Mw the W- 
boson mass, mt the top-quark mass and \i the 
renormalisation scale. S^m^/M^), tjb and 
Cs(/i) are perturbatively-calculated quantities. 

An alternative approach, in which many theo- 
retical uncertainties cancel, is to look at the ratio 
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is the meson mass. The matrix el- 
\ can be parameterised as 
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where fB d{s} is the decay constant, and -Bs d(s) the 
£?-parameter of -B°( s ) mesons. 

In this work, we obtain the ratio r s d from the 
direct calculation of Mbs/Mbd as well as from 
the calculations of /s s //s d and Bb s /Bb c1 - 

2. SIMULATION DETAILS 

We use the tadpole-improved Sheikholeslami- 
Wohlert (SW) quark action, 



_ qW 



E 



q{x)P llv {x)a IJ , v q{x) 



to perform simulations on a 24 3 x 48 lattice at 
/3 = 6.2 and a 16 3 x 48 lattice at /3 = 6.0. Here 
is the standard Wilson action, the bare 
gauge coupling, csw the clover coefficient, k the 
hopping parameter, and P M „ a lattice definition 
of the gauge-field strength tensor. Table [l] gives 
the simulation parameters. We use KLM normal- 
isation for the quark fields. 

3. OPERATOR MATCHING 

Matching onto the MS scheme is performed at 
one-loop in perturbation theory using the cou- 
pling cvss(/i) defined from the plaquette jjj . Since 
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Table 1 



Simulation parameters, kq and n q are the heavy- and light-quark hopping parameters. 





# configs. 


csw 






6.0 


498 


1.48 


0.13700 0.13810 0.13856 


0.114 0.118 0.122 0.126 0.130 


6.2 


188 


1.44 


0.13640 0.13710 0.13745 


0.120 0.123 0.126 0.129 0.132 



the clover-leaf interaction term is proportional to 
go, we can use the perturbative results obtained 
from a tree-level clover action Q with modifi- 
cations appropriate for tadpole-improvement and 
KLM normalisation. 

For the matching of four-fermion operators, we 
use the basis 



I a I 



la I 



Qlat 
Qlat 



7m x 7m +7^75 x 7^75, 
7m x 7a* - 7m 75 x 7^75, 

IX 1 + 75 X 75, 

/ x 7-75 x 75, 

O tnj X <7 nlJ . 



We set the coupling and matching scales to 
/i = — and, for consistency with the literature, 
run divergent operators to 5 GeV, using 2-loop 
continuum RG in the MS scheme with the appro- 
priate number of flavours. 

To estimate the systematic error associated 
with the one-loop matching, we vary the scale 
/j, in a range from 1/a to n/a. Decay constants 
are not affected since they are normalized by /„• 
and _B-parameters change by about 3% (f v varies 
by approximately 3%). Since we are mainly in- 
terested in SU (3)-breaking ratios for which these 
effects are even smaller, we neglect these small 
variations in what follows. 

4. ANALYSIS AND RESULTS 

We determine k c and k s from pseudoscalar me- 
son masses. We set the scale with M p for spec- 
tral quantities and f w for decay constants. In 
fact, these two quantities yield remarkably similar 
scales. (See Table |[) We then linearly extrapo- 
late and interpolate heavy-light decay constants, 
B-parameters and AB = 2 matrix elements to 
k c and k s , keeping k c , k s , aM p and af v in the 
bootstrap loop. Fig. 1 shows examples of these 
extrapolations. 



Table 2 

Critical and strange hopping parameters and in- 
verse lattice spacings. 





6.0 


6.2 


a~ 1 (M p )(GeV) 
a- 1 (/ w )(GeV) 


0.13924(1) 
0.13757(8) 
1.96(5) 
1.92(4) 


0.13793(1) 
0.13670(9) 
2.57(8) 
2.58(9) 



= 6.2 



0.120 
0.126 
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Figure 1. Light-quark-mass dependence of the 
heavy-light B-parameter, Bp, and extrapolation 
(interpolation) to m = k c (k s ) at /3 = 6.0 and 
6.2. 



For heavy-quark (HQ) extrapolations, we de- 
fine (Mp is the heavy-light meson mass) 



n P> Z A V M p \a s (M B ) 



*af= 2 (M p ) 



a^MMp \ a s {M P )\ 11 



M P {a s (M B ) 



Then for X(M P )=^ f (Mp),^ AF=2 (Mp), B(M P ) 
and 5f7(3)-breaking ratios, HQET predicts 

X(Mp) = A x {l + B x (^) + C xi ^r + . " 

Fig. 2 shows examples of the HQ extrapolations. 
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HQ Scaling 



□ :B — meson 
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Figure 2. HQ scaling of $/ and B-parameter 
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Figure 3. Chiral and HQ extrapolation of the 
matrix element. 



For 5'[/(3)-breaking ratios, we find that tak- 
ing the ratio before or after the HQ extrapo- 
lation leads to nearly indistinguishable results. 
We use the former for our final results since 
5?7(3)-breaking ratios have milder HQ-mass de- 
pendences. 

Our main results are summarised in Table [|. 
We obtain r s d from the direct calculation of 
M bs /M b d as well as from j Ba / fs d and B Bs / B Bd ■ 
Our results for the direct calculation are consis- 
tent with those of ||, obtained with propagating 
Wilson quarks, and, at f3 = 6.0, with the static 
result of j|. However, as Fig. 3 suggests, it is 
more difficult to control the chiral and HQ ex- 
trapolations of the matrix elements in the direct 
calculation because these extrapolations are more 
pronounced. 

Because we have results at only two values 
of the lattice spacing, we cannot extrapolate to 



Table 3 

Summary of results. Errors are statistical only. 
I 

Jd (MeV) 
/fl.(MeV) 
/ b (MeV) 

fDs 

it 

S^°(5GeV) 
B B fo (5GeV) 



Be 

MB* f Bs 



Me 



2 B Bs 
B b 



6.0 


6.2 


239(6) 


221(9) 


213(6) 


193(10) 


221(7) 


190(12) 


191(10) 


161(16) 


1.12(1) 


1.15(4) 


1.15(4) 


1.18(8) 


0.86(2) 


0.85(2) 


0.83(4) 


0.85(3) 


1.03(3) 


0.99(3) 


1.38(7) 


1.37(13) 


1.52(19) 


1.70(28) 



the continuum limit. We therefore consider the 
(3 = 6.2 results to be our best, noting that de- 
cay constants may still suffer from relatively large 
discretisation errors (roughly a 2a effect between 
6.0 and 6.2) while S'C/(3)-breaking ratios and B- 
parameters are consistent within errors at the two 
(3 values. 

Although formally one need not include the 
ad^P correction to the axial current when using 
a mean-field improved, tree-level clover action, it 
would be interesting to investigate its effect on 
our results in view of understanding how non- 
perturbatively, 0(a)-improved decay constants 
may behave. We plan to do so in the future. 

For a comparison of our results with other re- 
cent results, we refer the reader to 
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